Degeneracy measures for the algebraic classification of numerical spacetimes 
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We study the issue of algebraic classification of the Weyl curvature tensor, with a particular focus 
on numerical relativity simulations. The spacetimes of interest in this context, binary black hole 
mergers, and the ringdowns that follow them, present subtleties in that they are generically, strictly 
speaking, Type I, but in many regions approximately, in some sense, Type D. To provide meaning 
to any claims of "approximate" Petrov class, one must define a measure of degeneracy on the space 
of null rays at a point. We will investigate such a measure, used recently to argue that certain 
binary black hole merger simulations ring down to the Kerr geometry, after hanging up for some 
time in Petrov Type II. In particular, we argue that this hangup in Petrov Type II is an artefact 
of the particular measure being used, and that a geometrically better-motivated measure shows a 
black hole merger produced by our group settling directly to Petrov Type D. 

PACS numbers: 04.25. D-,04.20.Cv,04.20.Gz,04.25.dg 



I. INTRODUCTION: 

The marvelous improvements in the technology of nu- 
merical relativity in recent years present opportunities for 
revolutionizing our understanding of the classical gravi- 
tational field. In the past, much of this understanding 
has come from studying solutions with extreme symme- 
try, and perturbations of such solutions. However, with 
the help of numerical methods, truly generic simulations, 
particularly of multiple black hole systems, can now be 
carried out in full general relativity. 

While this work is undertaken, one must keep in mind 
the fundamental nature of general relativity and its so- 
lutions. In particular, the general covariance of the the- 
ory is not naturally reflected in the numerical context, 
where gauge fixing is fundamentally required in the form 
of coordinate and tetrad choices. In practice, such gauge 
choices are tailored to numerical convenience (or neces- 
sity), rather than to physical relevance. Such a simple 
task as checking that a black hole merger settles down to 
a Kerr geometry can be clouded by the arbitrariness of 
the simulation coordinates. 

One way of dealing with these ambiguities would be 
to apply coordinate transformations to numerical simu- 
lations a posteriori to represent these spacetimes in phys- 
ically preferable coordinates, if they exist. If one needs 
to map all quantities to an entirely new coordinate grid, 
then some accuracy would presumably be lost to the in- 
terpolation process, especially if changes of the time func- 
tion require interpolation in time. More important, how- 
ever, is the difficulty of fixing physically preferred coor- 
dinate systems in strongly dynamical and nonsymmetric 
spacetimes at all. 

Another, perhaps complimentary, approach is to fo- 
cus physical analysis on partially (or if possible, totally) 
gauge invariant quantities. For example, a major tool 



in the analysis (and construction) of exact solutions in 
general relativity is the algebraic classification system of 
Petrov and Pirani [ll-Q , in which the Weyl tensor at any 
given point in spacetime is classified according to the al- 
gebraic properties of its associated eigenbivector problem: 



C a b cd X cc [ — AX a b 



(1) 



Another view of this classification system, with a more 
geometrical flavor, was expounded particularly by Bel [3] 
and Penrose Q . In this approach one classifies the Weyl 
tensor in terms of the degeneracy of the so-called princi- 
pal null directions, null vectors defined up to scale by the 
equation: 



k e k f k[ a C b ] ef [ c k d ] = 0. 



(2) 



One can show (most easily in spinor language) that this 
equation is always satisfied by exactly four null rays, 
counting multiplicities. If all four of these null directions 
are distinct, the spacetime is said to be algebraically gen- 
eral or Type I at that point in spacetime. If two of them 
coincide, the spacetime is said to be Type II there. If 
three, Type III. If all four principal null directions co- 
incide, the spacetime is said to be Type N, or null, in 
analogy with the pure radiation fields of vacuum electro- 
dynamics. If the principal null directions coincide in two 
distinct pairs, then the spacetime is said to be Type D. 
The Kerr and Schwarzschild geometries are famous ex- 
amples of globally Type D spacetimes, so in some sense, 
one may hope to infer that a spacetime is "settling down 
to an approximate Kerr geometry" if its Petrov type "set- 
tles down" to Type D (assuming that one has ruled out 
other, non-Kcrr, Type D spacetimes). 

This line of reasoning was taken up in a paper by Cam- 
panelli, Lousto, and Zlochower [(|. The central tool in 
their approach is a certain complex polynomial equation: 



* 4 A 4 + 4* 3 A J + 6* 2 A 2 + 4*iA + * = 0, 



(3) 
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the degeneracy of whose roots is known to correspond to 
the degeneracy of the principal null directions (assuming 
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that ^4 is nonzero). The coefficients ^ of this polyno- 
mial are the so-called Weyl scalars, components (defined 
in Eq. (j9|) below) of the Weyl tensor in a given Newman- 
Penrose null tetrad. If ^4 is nonzero, then the funda- 
mental theorem of algebra ensures that the polynomial 
has exactly four complex roots, counting multiplicities. 

Once the four roots have been computed for Eq. ^ 
at any given point, then one can also compute six distinct 
positive-definite root differences: 

Afc-HAi-Ajl. (4) 

If two of these root differences vanish and the other 
four are nonzero, meaning that the four roots coincide 
in two distinct pairs, then the spacetime is Type D at 
that point. In Rcf. [6[ Campanelli et al. took the next 
logical step: interpreting the two smallest Aj 3 - values as 
measures of the "nearness" of an algebraically general 
(in that case, numerical) spacetime to Petrov Type D. 
While this is a reasonable interpretation of Ay and we 
will not suggest any fundamental modification to this ap- 
proach of defining approximate algebraic speciality, there 
are important subtleties in this interpretation, not fully 
explored in Ref. @. These subtleties relate to the ge- 
ometrical meaning of Ay and its behavior under tetrad 
transformations. The main purpose of this paper is to ex- 
plore these subtleties, present an alternative degeneracy 
measure that avoids certain blowups that arc intricately 
related to the choice of tetrad (and should therefore not 
be considered physically relevant), and apply both de- 
generacy measures to a numerical simulation from the 
SpEC code @- In the process we will also investigate an 
interesting conclusion from Ref. [(| : that in the ringdown 
of a binary black hole merger to Kerr, the spacetime ap- 
proaches Petrov Type II very quickly, and Type D much 
later. We will argue that this conclusion is due essentially 
to a coordinate singularity on the space of null rays, and 
the fact that the tetrad used in Ref. @ was much bet- 
ter suited to representing the degeneracy of one pair of 
principal null directions than the other pair, when the de- 
generacy measure Ay is used. The alternate degeneracy 
measure that we will introduce, 0y defined in Eq. (|45[) 
below, shows both pairs of principal null directions ap- 
proaching degeneracy at the same rate. 

Though much of the discussion in this paper centers 
upon the behavior of these measures of degeneracy under 
tetrad transformations, we will unfortunately not be able 
to provide a measure of nearness to Petrov Type D that 
is fundamentally any more invariant than Ay. This is 
because no such measure appears to exist. Geometrically 
this fact can be understood in terms of the nonexistence 
of a boost-invariant geometry on the space of null rays 
in Minkowski space, an issue referred to physically as the 
"relativistic aberration of starlight." This viewpoint is 
explored in more detail in Section [IVI below. 

The issue can also be understood at the algebraic level, 
as in Petrov's original construction. The problem shown 
in Eq. ([T]) can be written more compactly if one works 
in the three-complex-dimcnsional space of anti-self-dual 



bivectors rather than in the six-real-dimensional space of 
real bivectors. In this space, the eigenbivector problem 
can be written as: 

W ab cd Z cd = AZ ab , (5) 

where W abcd := C abcd + i *C abc d, *Z ab = —iZ ab , and A is 
a complex number. Because this is a three-dimensional 
problem one can expect three possible values for A, 
though the fact that the Weyl tensor is tracefree implies 
that these three eigenvalues must sum to zero. The de- 
generacy of the eigenvalues and the completeness of the 
corresponding eigenspaces determine the classification of 
the Weyl tensor at the point under consideration. If all 
three eigenvalues are distinct then the spacetime is alge- 
braically general. If two roots coincide, then the space- 
time is either Type II or Type D. If all three coincide 
(and therefore vanish, as they must sum to zero) then 
the spacetime is either Type III, Type N, or conformally 
flat. The eigenvalues are geometrically defined at each 
point in spacetime, independent of the vector basis used 
to represent the eigenproblem. The differences between 
these eigenvalues can therefore be used to construct in- 
variant measures of the approach to algebraic speciality. 
For example, the absolute value of the difference between 
the two nearest eigenvalues can be thought of as such an 
invariant measure. Unfortunately this measure isn't very 
specific: it vanishes for Petrov Types II, D, III, and N. 1 
The latter pair can be distinguished from the former pair 
by the fact that all three eigenvalues vanish in Type III 
and Type N, but distinguishing Type II from Type D, or 
Type III from Type N, requires more information than 
just the eigenvalues. 

If two of the eigenvalues in Eq. ([3]) coincide, so that 
the eigenvalues can be written as {A, A, — 2A}, then the 
distinction between Petrov Type II and Type D can be 
made by the following quantity 2 : 

T ab cd := (W ab ef - Alf f )(W ef cd + 2AI e J), (6) 

where I^ d is the identity operator on the space of anti- 
self-dual bivectors. The object T ab cd vanishes in Type 
D, but not in Type II The difficulty with using this 
as a measure of nearness to Petrov Type D is that it is 
a tensorial object, and its components are, by definition, 
basis-dependent. In order to collapse this object to a 
single number for each point in spacetime, one might 



1 In this sense it is like other scalar measures of algebraic special- 
ity, such as the "cross ratio" of principal null directions, defined 
in Ref. whose explicit relationship to the eigenvalues is de- 
scribed in Sec. 8.3 of Ref. or the Baker-Campanelli "speciality 
index" [Toll , which takes a special value for any type of algebraic 
speciality, but cannot distinguish between the various types. 

2 Incidentally, if one wishes to avoid the assumption that the space- 
time is at least Type II, such that the eigenvalues can be written 
as {A, A, — 2A}, this can be done with the help of certain curva- 
ture invariants. See Ref. [Till . 
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hope to construct a positive-definite tensor norm: 



II. TETRAD DEPENDENCE 



Q := m ae m bf m c a m dh T ab cd T &S gh , (7) 



where m ab is a positive-definite inner product on space- 
time. Unfortunately, the only inner product that one nat- 
urally has available on spacetime is the indefinite space- 
time metric. If a timelike "observer" is introduced, with 
unit tangent vector u a , then one can construct a positive 
definite inner product as: 



m a b ■— gab + 2u a Ub, 



(8) 



but then the quantity Q is not strictly a scalar, as its 
definition is dependent on the extra structure of this ob- 
server. 

Though the language is very different in the geometric 
approach involving principal null directions, we will find 
in Section IIVI that the ambiguity in defining a measure 
of "nearness" to Petrov Type D is in that context essen- 
tially the same as here, requiring the choice of a timclikc 
observer at every point in spacetime. While this state 
of affairs seems to endanger any attempt at defining the 
nearness to any specific Petrov class, there are some cases 
where a well-defined fleet of observers can be chosen. In 
particular, in any stationary spacetime, one can choose 
the stationary observers. In cases such as the ringdown 
to Kerr geometry, one can expect an "approximate" sta- 
tionarity to be approached at late times, again provid- 
ing a preferred class of observers at least during the late 
ringdown. A major practical goal of this paper will be to 
study this ringdown process, as in Ref. [6[. In particu- 
lar, we will argue that the degeneracy measure Ay that 
was used in Ref. @ is in some sense adapted to a null 
observer that happened in that case to be nearly aligned 
with one of the nearly degenerate pairs of principal null 
directions, making this pair of null directions seem much 
more degenerate, and the other, much less. This causes 
the appearance of a holdup in Petrov Type II before the 
spacetime geometry falls to Type D. 

The structure of this paper is as follows: in Sec. HT1 we 
will investigate the ambiguity of the measure Ay under 
tetrad rotations, particularly those that leave the time- 
like tetrad leg fixed. In Sec. IIIII we will emphasize the 
fact that a tetrad well-suited to gravitational wave ex- 
traction, in particular the quasi-Kinnersley tetrad (l2T |. 
may be particularly ill-suited to measuring the nearness 
to Petrov Type D using Ay. In Sec. IIVI we will describe 
the geometry underlying Ay in spinorial language, and 
in the process motivate a modification that is much bet- 
ter suited to situations such as the ringdown to Kerr 
geometry. In Sec. [V] we will present numerical results 
applying these degeneracy measures to a binary black 
hole merger simulation, demonstrating in detail the ap- 
proach to Petrov Type D. Finally in Sec. IVII we conclude 
with further discussion of the subtleties that have been 
addressed, and those that remain. 



The method put forth in Ref. |6| to define nearness 
to a Petrov class begins with the polynomial in Eq. ([3]), 
whose coefficients are components of the Weyl tensor in 
a Newman-Pcnrose tetrad UM : 



*o := C abcd £ a m b £ c m d , 
*i := C abcd tn b £ c m d , 



(9a) 
(9b) 



* 2 := -C abcd {l a n h l c n d -l a 

*3 := C abcd n a £ b n c m d , 
*4 := C abcd n a m b n c m d . 



n b m c m d ) , (9c) 



(9d) 
(9e) 

The tetrad {£ a ,n a ,m a ,rri a } is made up of two future- 
directed real null vectors £ a and n a and two complex 
conjugate null vectors to° and rn a with spacelike real 
and imaginary parts. These vectors are normalized by 
the conditions: 



£ a n a = -1, 
m a m a = 1, 
£ n m a = n n m a = 0. 



(10) 

(11) 
(12) 



These normalization conditions are preserved by three 
types of tetrad transformations which, taken together, 
are equivalent to the proper Lorentz group. First, there 
are the "null rotations about £ a " sometimes referred to 
as the "Type I" transformations 3 : 



to i— > to 
n a h-> n a + am a + arrf + oat, 



(13a) 
(13b) 
(13c) 
(13d) 



where a is a complex number, and can vary over space- 
time. Second, there are the null rotations about n a , 
sometimes referred to as "Type II" transformations: 

T i-> £ a + bm a + bm a + bbn a , (14a) 

m a h-> m a + bn a , (14b) 

m a ^ m a + bn a : (14c) 

n a ^ n a , (14d) 

for complex b. Third, there are the "spin-boost" transfor- 
mations, sometimes referred to as the "Type III" trans- 
formations: 



i — ^ c 



2na 



m a ^ e 2iarg( C ) m a 5 
—a ^ e _2iarg(c)— 

n a i y \c\- 2 n a , 



(15a) 
(15b) 
(15c) 
(15d) 



3 To avoid confusion with the Petrov types, we will hereafter re- 
fer to tetrad transformations as "null rotations about £ a ," "null 
rotations about n a ." or "spin boosts." rather than "Type I," 
"Type II," or "Type III." 
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for complex c. 

These transformation laws for the tetrad imply trans- 
formation laws for the Wcyl scalars. Under the null ro- 
tations about £ a , Eqs. (fl~3)). the Weyl scalars transform 
as: 



*i i y * 1 + a* , 

* 2 ^ *2 + 2a*i 

* 3 i y * 3 + 3a^ 2 + 3a 2 *i 

* 4 i y * 4 + 4a^ 3 



6 a 2 * 2 



a 3 * , 
4a 3 *i 



(16a) 
(16b) 
(16c) 
(16d) 



-a 4 * - (16e) 



Under null rotations about n a , Eqs. (HU), the Weyl scalars 
transform as: 



*o !->■ & *4 



46 3 *, 



66 2 * 2 +4fe*i + *o, (17a) 



*i i y 6 3 * 4 + 36 2 * 3 + 36* 2 

* 2 & 2 *4 + 26*3 + *2, 

* 3 I y 6*4 + * 3 , 

*4 h-» * 4 . 



*1, 



(17b) 
(17c) 
(17d) 
(17e) 



Finally, under the spin boosts, Eqs. (|T5| . the Weyl scalars 
simply rescale, as: 



*„ h-> c 2(2 - n) *„. 



(18) 



The transformation laws for the coefficients of the poly- 
nomial in Eq. (|3|) imply transformation laws for the roots. 
It is straightforward to show that under the transforma- 
tion in Eq. (|16p . the roots of the polynomial transform 
as: 



A i-» 



A 



aA + 1 



(19) 



Under transformations of the form ([TT]) , the roots trans- 
form as: 



A h- > A + 6. 



(20) 



Finally, under spin-boost transformations, Eq. (|18|) . the 
roots transform as: 



A H. c 2 A. 



(21) 



In Ref. [6j, nearness to Petrov Type D was mainly ar- 
gued through the approach of the absolute values of root 
differences (Ay as defined in Eq. (0|) to zero. While this 
quantity would indeed be expected to vanish when Xi 
and Xj constitute a degenerate root pair, if they are not 
exactly degenerate, then the foregoing discussion implies 
that this difference is not invariant under tetrad transfor- 
mations. The transformation in Eq. (|20p would leave Ay 
unchanged, but that in Eq. (f2"Tj) would directly rescale 
any given root difference (though the complex phase of c 
would not appear in the absolute value) , and transforma- 
tions of the form (|19[) would change Ay in a more com- 
plicated way. Arbitrary Lorentz transformations, given 



(£ a + n a 
V2 Re [m a ] , 
V2 Im [m a ] , 
(£ a - n a )/V2. 



by arbitrary combinations of the above transformations, 
could alter |A, — Xj\ in a very complicated manner. 

To investigate the practical relevance of this tetrad 
ambiguity in the degeneracy measure Ay, let us con- 
sider a particular case of possible physical relevance that 
requires a combination of all three of the above tetrad 
transformations. Take the case where one has a particu- 
lar timclikc vector defined at a point in spacetime, for ex- 
ample a timelike Killing vector, or a kind of approximate 
Killing vector generating time translations in a spacetime 
that is approaching stationarity in some sense. Given a 
Newman-Penrose tetrad {£ a , n a , m a , fn a }, one can con- 
struct a standard orthonormal tetrad in the following 
way: 

(22a) 

-l •— v ^ i« [in j , (22b) 
e a 2 := V2 Im [m a ] , (22c) 

(22d) 

Rotations of the tetrad legs in the e*i-e2 plane are easily 
accomplished, through a simple spin-boost transforma- 
tion with the parameter c = e 4 */ 2 . Such rotations also 
however leave the degeneracy measure Ay unchanged. 
For a nontrivial test case, consider a rotation in the e\- 
(=3 plane: 

e ' = e (23a) 
ei = cos($)ei - sin($)e 3 (23b) 
e 2 ' = e 2 (23c) 
e 3 ' = cos($)e 3 + sin($)e*i (23d) 

A straightforward calculation shows that such a trans- 
formation can be carried out by a sequence of the above 
transformations. First, one makes a null rotation about 
£ a , Eqs. (fT3]) , with parameter a = — tan($/2). Second, 
there is a null rotation about n a , Eqs. (fT4")l , with param- 
eter 6 = (1/2) sm($). The final step is a spin boost, 
Eqs. (|15p . with parameter c = sec($/2). In this partic- 
ular case, all three parameters are real. Composing the 
transformation laws for the roots, Eqs. (|19p - (|2"Tj) , with 
these parameters, the resulting transformation law is: 



A' 



Acos($/2) - sin($/2) 
Asin($/2) +cos($/2)' 



(24) 



If we express A as a ratio of two complex numbers, A 
£/?7, then Eq. (|24p takes a very simple matrix form: 



cos($/2) -sin($/2) 
sin($/2) cos($/2) 



(25) 



The general form of this matrix, for arbitrary reorienta- 
tions using three Euler angles, is given in Eq. (1.2.34) of 
Ref. Q. The SL(2,C) form of this transformation sug- 
gests a spinorial interpretation of A, a point to which we 
will return in Sec. IIVI 
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tetrad rotation angle <I> 

FIG. 1: Behavior of the degeneracy measure A„ under the 
tetrad rotation in Eq. (|23l) for a particular (though essentially 
arbitrary) choice of roots, stated in the text. Under a rota- 
tion through 180 degrees, the root pair that originally seemed 
more degenerate becomes less degenerate, and the pair that 
originally seemed less degenerate becomes more degenerate. 



For now let us consider the behavior of the degeneracy 
measure Ajj under these spatial rotations. For concrete- 
ness, consider the case where the four roots of Eq. ([3]) are 
A, = .005 + .047*, A 2 = .005 + .05i, A 3 = -5 + 15i, and 
A4 = — 5 + 15.5i. These values are chosen to very roughly 
mirror the late-term values seen in Fig. 8 of Ref. Q, with 
degeneracies roughly similar to those seen in Figs. 3 and 
4 of that paper. The values estimated here are extremely 
rough, and should not be taken as having any quantita- 
tive importance, but merely as tools for illustrating the 
qualitative features of the transformation law in Eq. (f24)) . 
So long as one pair of nearly-degenerate roots is larger, 
by a few orders of magnitude, than the other pair, the 
qualitative behavior that we will describe seems roughly 
the same regardless of the particular choice of roots. 

The degeneracy measure Ay for the two most nearly 
degenerate root pairs, under rotation of the ei-e^ plane, 
is shown in Fig. [T] If the tetrad's spatial legs were ro- 
tated through about ninety degrees, then both root pairs 
would appear equally close to degeneracy. If the tetrad 
were rotated through 180 degrees, then the root pair that 
originally appeared closer to degeneracy would begin to 
seem farther away from it, and the one that originally 
seemed less degenerate would seem more so. 

This variation in the degeneracy measure can be in- 
terpreted as a coordinate effect. The quantity A has no 
inherent geometrical meaning without a particular refer- 
ence tetrad. It is essentially a coordinate on the space 
of null rays at a point in spacetime. This space of null 
rays is topologically a two-dimensional sphere, as can be 
demonstrated by cutting a future null cone with a space- 
like hypcrplane, such as the t = 1 plane in Minkowski 



space. A two-sphere cannot be covered smoothly with a 
single coordinate patch. If the quantity A is taken as a 
(complex) coordinate labeling all the null rays at a point, 
then there must be a coordinate singularity somewhere, 
near which coordinate distances are particularly ill-suited 
to representing the true geometry that may be defined on 
the manifold. We will study this issue in more detail in 
Sec. IIV1 For now we simply note that the locations of 
the sharp peaks in Fig. [T]seem to imply that such a co- 
ordinate singularity may have a particularly strong effect 
in the original, unrotated tetrad. The following section 
gives an extreme example of this effect. 

III. THE QUASI-KINNERSLEY TETRAD: 

It appears from the results of the previous section that 
a tetrad that seems reasonable for purposes of wave ex- 
traction can be particularly ill-suited to the problem of 
defining nearness to a Petrov class. To investigate this 
point in more detail, here we consider a special family of 
tetrads designed especially for wave extraction. 

Consider an algebraically general spacetime (eventu- 
ally we will allow this spacetime to "asymptote" toward 
Petrov Type D, but we will consider it always to be, 
strictly speaking, Type I). As described in Ref. at 
any point where the Weyl tensor is Type I, there are 
precisely three distinct families of tetrads in which two 
particular Weyl scalars vanish, Vti = ^3 = (they each 
amount to families of tetrads, rather than three partic- 
ular tetrads, because this condition is preserved by the 
spin-boost freedom). A particular tetrad field, chosen 
from these three families to coincide with the conven- 
tional Kinnersley tetrad near infinity, is often referred to 
as a quasi- Kinnersley tetrad. The usual purpose of such 
a tetrad is to aid in gravitational wave extraction, where 
the relative uniqueness of the tetrad provides a preferred 
reference frame in which to define gravitational radiation. 
Such a tetrad also simplifies the polynomial in Eq. 

* 4 A 4 +6* 2 A 2 + * = 0. (26) 

If, as we are assuming, the spacetime is strictly Type I, 
and not a more special algebraic type, then "J 4 and v&o 
will be nonzero. In the limit that the spacetime asymp- 
totes to Type D, they will both settle to zero, indicating 
a failure of the polynomial roots to represent the princi- 
pal null directions in the conventional sense. What we 
wish to investigate is the behavior of these roots as this 
limit is approached. 

Carrying on under the assumption that $4 is nonzero, 
the roots of Eq. (|26|) arc readily found. 

If we now consider the approach to a Kerr geometry, 
in which the quantity ^o^/C^l) approaches zero, we 
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can expand the square root in the above expression to 
first order in this small quantity 4 : 



A 2 




(28) 
(29) 



So in the Kerr limit, as \&4 — > and "to — > 0, two of these 
roots approach zero, and so does their difference, but the 
other two approach infinity (this is a standard behavior 
of polynomial roots as the leading polynomial coefficient 
approaches zero). Moreover, they approach the point at 
infinity from different directions, so their difference also 
approaches infinity. Geometrically, one would think that 
the problem is solved if the roots are considered not as 
numbers on the complex plane, but as points on the Rie- 
mann sphere. The roots that blow up would then be 
taken as approaching a degenerate root at the point at 
infinity. In the following section, we will motivate such 
a viewpoint in detail, and in the process, outline the ge- 
ometrical meaning of the degeneracy measure Ajj and 
present an alternative that avoids the danger of repre- 
senting any particular null ray as a "point at infinity." 

Before moving on, though, we should investigate 
the robustness of this behavior under tetrad rota- 
tions. In practice, the tetrads used in numerical 
relativity simulations arc usually simple coordinate- 
adapted tetrads, rather than carefully constructed quasi- 
Kinnersley tetrads. But because they are usually adapted 
to a spherical coordinate basis, they very roughly tend 
to approximate the quasi-Kinnersley tetrad during black 
hole ringdown, by force of topology alone. For this rea- 
son, it is interesting to investigate the behavior of the 
polynomial roots not only in the quasi-Kinnersley tetrad, 
but also in tetrads slightly offset from it. 

In particular, consider the ringdown to a Kerr black 
hole, where in the true Kinncrsley tetrad of a Kerr back- 
ground one would expect the absolute value of $4 to 
approach zero exponentially in time at a rate determined 
by the quasinormal frequencies of the hole. The roots 
± \f— 6^2/^4 would then be expected to grow exponen- 
tially at half that rate. Consider, for example, a case 
where y/— 6^2/^ a = icxp(r) for some time function r. 5 
The roots ±iexp(r), if the tetrad were rotated spatially 
as in Eq. (J24J) , would instead take the values: 



A'.. 



±zexp(r) cos($/2) — sin($/2) 
±i exp(r) sin($/2) + cos($/2) ' 



(30) 



Note that the numerator in this quantity, ^o^i which we are 
evaluating in a "transverse frame" — one where = ^3 = 
— is the Beetle-Burko "radiation scalar" described in Ref. [13 . 
The imaginary factor i is inserted to avoid the rotated tetrad 
vector exactly coinciding with a principal null direction at some 
time, a possibility that, however possible, would not be expected 
to occur generically. 




-1 1 

tetrad rotation angle <I> 



FIG. 2: Profiles of the behavior of the degeneracy measure 
Aij under tetrad rotations of the form in Eqs. (I23|) from a 
quasi-Kinnersley tetrad, for various values of a fiducial time 
coordinate, assuming that this degeneracy measure grows ex- 
ponentially in this fiducial time coordinate for the true quasi- 
Kinnersley tetrad ($ = 0). The peak value grows exponen- 
tially in time, by construction, but values well outside the 
peak decay exponentially in time. The peak sharpens as it 
grows, so that values slightly offset from the peak grow ini- 
tially, and decay later. 



So in the limit that r — > 00, these roots would become 
degenerate at the value cot($/2), and their difference, 
as measured by A^ , would eventually fall to zero. The 
details of how this occurs are plotted in Figs. [2] and [3] 

In Fig. [2 the profile of A', := |A' + — X'_ |, as a function 

of tetrad rotation angle $ in Eq. (|5U|) . is shown for a 
few values of the fiducial time label r. Each curve is 
peaked, as in Fig. [TJ at the quasi-Kinnersley tetrad. The 
value of this peak grows exponentially in r, while the 
values well outside the peak (representing more arbitrary 
tetrads) decay exponentially in r. What is of particular 
interest to us is the behavior near $ = 0. Because the 

peak sharpens as it grows, values of A', slightly offset 

from $ = grow initially, and eventually decay. This 

behavior is more clearly visible in Fig. [3J where A^ 

is shown as a function of r for various choices of the 
offset angle. For any fixed nonzero value of $, the curve 
initially grows exponentially before eventually falling at 
the same rate. The smaller this rotation angle is, the 
later the curve turns around. So the nearer a tetrad 

is to quasi-Kinnersley, the longer it takes for A^ to 

eventually decay as one might naively expect. 

Incidentally, we should note that the "peak" in Fig. [5J 
and indeed in Fig. [TJ is actually a saddle point if con- 
sidered in a larger space of tetrad rotations. To keep 
the discussion simple, in Sec. [II] we considered only rota- 
tions in the local e*i-e3 tangent plane. Had we considered 
the case of general rotations of the spatial tetrad, as in 
Eq. (1.2.34) of Ref. @, we would have found that the 
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FIG. 3: Behavior over time of the degeneracy measure 
AV_($) := |A' + - A'_|, for roots A' ± given by Eq. (O, for 
a few values of the tetrad rotation parameter $. Each curve 
initially grows exponentially in the time parameter r before 
eventually falling. The more offset the tetrad is from the 
quasi- Kinnersley tetrad (4> = 0), the sooner this turnaround 
occurs. 



degeneracy measure Ajj becomes infinite whenever the 



rotated tetrad 
direction. 



vector coincides with a principal null 



IV. INTERPRETING DEGENERACY 
MEASURES: 

The geometrical underpinnings of the polynomial in 
Eq. and the sense in which A constitutes a coordi- 
nate on the space of null rays, are most cleanly explained 
in the language of two-component spinors. Because many 
numerical relativists are unfamiliar with this formalism 
I will attempt to keep the discussion self-contained by 
briefly reviewing crucial elements as we go along. For a 
detailed account of spinor methods in spacetime geom- 
etry, see Refs. @, @, or for a more compact treatment 



specifically geared to numerical relativists, see Ref. [15 1. 

Throughout this paper, objects with capital latin in- 
dices will be referred to as spinors, elements of a two- 
complex-dimensional vector space (or its higher tensorial 
orders). The complex conjugate of a spinor is also a 
spinor, but is defined in a different spinor space, because 
complex conjugation does not commute with multiplica- 
tion by a complex scalar. To distinguish objects in spinor 
space from objects in the complex conjugate space, we 
will apply the standard convention of appending indices 
referring to the latter space with a prime: 



a A = a A ' . 



(31) 



pie correspondence exists between spinor space and 
Minkowski space (and therefore also to the tangent space 
to spacetime at any given point, given an orthonormal 
tetrad). From a spinor a A a unique vector can be con- 
structed in Minkowski space: 



jra „A—A'a 

V = a a a aa 1 



(32) 



where the a a AA' are soldering forms, specifically referred 
to as Infeld~van den Waerden symbols, conventionally 
represented as Pauli matrices. In practice, the transfor- 
mation provided by o~ a aa 1 Is often (and hereafter) taken 
as implied, with pairs of capital latin indices (one primed 
and one unprimed, with the same letter) taken to corre- 
spond abstractly to a single spacetime index. 

Vectors defined directly from univalent spinors as in 
Eq. (|3"2"|) turn out always to be null. For that reason uni- 
valent spinors can be understood as defining null vectors 
in spacetime. The standard geometrical interpretation of 
a univalent spinor (again, see Ref. [H, [l5| ) is as a "null 
flag," a null vector with a particular spacelike half-plane 
attached to it. This flag plane, encoded in the spinor's 
complex phase, is unimportant for our current purposes. 

The spacetime Weyl tensor can be written in terms of 
a four-index totally symmetric object called the Weyl 
spinor ^abcd and the antisymmetric metric cab on 
spinor space: 



W a 



bed 



* ABCD^A'B'^C'D' 



(33) 



Here, as in the introduction, W a bcd refers to the Weyl ten- 
sor in its complex anti-self-dual form, W a bcd '■= C a bcd + 
i * 'C a bcd- A basic result in spinor algebra (due essentially 
to the fundamental theorem of algebra) is that any to- 
tally symmetric spin tensor can be decomposed into a 
symmetrized product of univalent spinors. In particular, 
for the Weyl spinor, 



ABCD 



= o-^aPbIcSd)-, 



(34) 



for univalent spinors a a, Ab, 1c i $d defined up to arbi- 
trary complex scaling (any one of them can be scaled at 
the cost of inversely scaling another) . These are referred 
to as principal spinors of ^ abcd- Because the principal 
spinors are defined only up to a complex scaling, their 
corresponding null vectors arc defined only up to an ar- 
bitrary real scaling, and their flag planes are completely 
undefined. The corresponding null vectors, defined up to 
scale, are the principal null directions of the Weyl tensor 
at the spacetime point under consideration. 

Because the metric on spinor space, cab, is antisym- 
metric, all spinors have vanishing norm. 



OLAOi 



0. 



(35) 



For this reason, the condition for a A to be a principal 
spinor of the Weyl spinor is: 



Spinors are useful in relativity theory because a sim- 



^ABCDa A a B a c a D 



0. 



(36) 
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To consider this equation more concretely, we introduce 
a basis, {o A , l a }, in spinor space, normalized by the stan- 
dard condition cabo a l b = 1. Such a spin dyad is equiva- 
lent 6 to a Newman-Penrose tetrad through the definitions 
£ a = o A o A ' , n a = l a T a ' , m a = o A T A ' . Given a spin dyad, 
an arbitrary spin vector can be written as: 



a 



rjo 



(37) 



for complex components 77, £. Because we are only inter- 
ested in spinors up to arbitrary complex scaling, we can 
divide by 77 to let: 



(38) 



where £ = is a possibly-infinite complex number, 
an element of the onc-point-compactified complex plane, 
C U {00}, the Ricmann sphere. 

Scaling an arbitrary spinor to be of this form, and in- 
serting it into Eq. f|36[) . the resulting equation is: 



tf 4 <T + 4^ 3 C + 6* 2 C + 4*iC + *o = 0, 



(39) 



where we have used the standard spinorial definition of 
the Weyl scalars: 



*0 


= ^ABCDO A B C D , 


(40a) 


#1 


= ^ABCDO A B C L D , 


(40b) 


^2 


= * ABCD A B L C L D , 


(40c) 


#3 


= *ABCDO A L B l C L D , 


(40d) 


^4 


= *ABCDL A L B L C L D . 


(40e) 



We thus find, comparing Eq. ((3]) with Eq. (|3"9"|) . that 
the quantity A can be interpreted as the complex stereo- 
graphic coordinate £ on the Ricmann sphere, and in par- 
ticular, as defining a spinor a A of the form in Eq. (|38[) in a 
given spin dyad. Hereafter we will consider £ and A to be 
the same quantity, and use the symbols interchangably. 

This stereographic interpretation of £ (or A) is not 
merely a formality. As described in Chapter 1 of Ref. Q , 
the space of future-directed null rays at a point in space- 
time is topologically a two-sphere. This can be demon- 
strated by cutting a future null cone with a spacelike 
3-plane, given by t = 1 in some local Minkowski coordi- 
nate system. Furthermore, if we choose a particular such 
cut, whose intersection with the null cone we will label 
S + and call the anti- celestial sphere, after Ref. Q, the 
metric induced on this two-sphere from that in the local 
Minkowski spacetime is: 



ds 2 = 



(1 + CC) 2 ' 



(41) 



Strictly speaking, the correspondence is two-to-one, as the spin 
dyad {— o A , — l a } defines the same tetrad as {o A ,i A }. The dis- 
tinction, however, is not important here. 



where for the coordinate we have chosen the ( value of 
the spinor, of the form in Eq. (f3"5|) . whose associated null 
direction intersects S + . Applying the transformation to 
conventional spherical coordinates, 

C = cot (6»/2), (42) 

we arrive at the standard form of the unit sphere metric: 

ds 2 = dd 2 +sm 2 (d)d(j> 2 . (43) 

As a geometrical method for defining the nearness of 
two null directions to degeneracy, one can consider the 
metric (|43| on the anti-celestial sphere. If Q and Q 
are two roots of Eq. ([3]), then one can translate them 
to spherical coordinates (8i,4>i), {6j,<f)j) by inverting 
Eq. (|42[) . and then use the metric distance function on 
the unit sphere, given by the haversine formula as: 



Qij :=2 arcsin { (sin 2 [(0 t - 8j) /2] 



sin ttj sin t/o sin 



in 2 [(&-&)/2]) 1/2 } 



(44) 



This can also be written directly in terms of the stereo- 
graphic coordinates as: 



&ij = 2 arcsin 



10-01 



(1 + OOX1 + GO-) 



(45) 



As one can verify by a direct substitution of Eq. (|24|) , this 
degeneracy measure is invariant under spatial rotations 
of the form in Eq. (|23[) , or indeed any tetrad rotation 
that leaves the timelike tetrad leg invariant. 

We must stress, however, that even this is not a totally 
invariant measure of degeneracy. In fact, there are fun- 
damentally as many degrees of ambiguity in this measure 
as there are in \Q — Q\. The ambiguity in is encoded 
in the choice of cut one makes to the null cone in order to 
construct S + . This can be interpreted physically as a re- 
sult of the special relativistic effect known as "relativistic 
aberration of starlight," by which the inferred geometry 
of the celestial (or as in this case, anti-celestial) sphere is 
conformally mapped under Lorentz boosts. 

A geometrical interpretation of the degeneracy mea- 
sure Ay := \Q — 01 can l> e found in spinor space. Two 
spinors a A and a A are proportional — and therefore their 
associated real null vectors are proportional — if and only 
if their antisymmetrized product vanishes: 



[0:1,02] 



A B 



0. 



(46) 



It is tempting to use this quantity as a measure of the de- 
generacy of the null rays associated with a A and a A , but 
we must remember to account for the scaling ambiguity 
of the spinors. If [0:1,0:2] is nonzero, then an arbitrary 
rescaling of cither spinor, which should not alter any rea- 
sonable measure of the degeneracy of the null rays, would 
directly rescalc [01,02]. This ambiguity must be fixed 
by imposing a condition on the scaling of a A and a A . 
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One possibility, given a particular spin dyad {o A , t A }, 
is to assume that the spinors are of the form (|38[) . with 
af = o A + di A , and a A = o A + ( 2 t A ■ This condition 
can be stated for the associated null vectors and V 2 a 
in terms of the Newman-Penrose tetrad as: 

V> a = -1, (47) 

for i e {1,2}, along with the conditions that the V? 
are real null vectors, and where the Newman-Penrose 
tetrad vector n a is defined from the dyad spinor i by 
n° := l a T a . This subset of the future null cone can be 
visualized as its intersection with a null hyperplane de- 
fined by t + z = \/2 in the local Minkowski coordinates. 
To the mind accustomed to Euclidean geometry, this in- 
tersection might be assumed to be a paraboloid. How- 
ever, interestingly, the Lorentzian structure of the space- 
time metric causes the intersection to be, in terms of the 
induced metric, a fiat two-dimensional plane, with Q a 
standard complex coordinate on this plane. In fact, the 
absolute value of the degeneracy measure [a±, a 2 ], under 
this particular normalization condition, is precisely the 
quantity A12 = | Cl_^ C2 1 • For this reason, the degeneracy 
measure used in [6| can be understood as a geometric 
distance between the two associated principal null direc- 
tions along the cut made by a null hyperplane through 
the future null cone. 

The degeneracy measure 8^ introduced in Eq. (|4"5)l 
can similarly be understood in terms of the symplectic 
product [01,02]. If the condition on the null vectors as- 
sociated with the a A is taken to be that Vi ■ §0 = — 1, for 
a timelike tetrad vector defined from a Newman-Penrose 
tetrad through Eqs. rather than V% ■ n = —1, then 

the spinors must be scaled as: 

af = ■ 1 {o A + Qi A ). (48) 

In this case the absolute value of [01,02] becomes: 

|[o 1; o 2 ]j = ^= IC ^~ C21 _ , (49) 

^/(1 + CiCi) (1 + C2C2) 

essentially equivalent to 0i2, as defined in Eq. (|4"5)l . 

The distinction between the degeneracy measures Ay 
and Qtj can therefore be understood as a distinction be- 
tween two different realizations of the geometry of the 
space of null rays at a point. If the geometry in this 
space is inferred by cutting the null cone with a null hy- 
perplane, the distance function on the set of null rays is 
given by Ay. If the cut is taken by a spacelike hyper- 
plane, the distance is given by 0,j. 

The ambiguity of these distance functions stems from 
the ambiguity of these cuts. Fundamentally there are 
equally many degrees of ambiguity in both types of cut. 
A spacelike hyperplane can be boosted in any of three 
directions, translating into three continuous degrees of 
ambiguity for Oy at each point in spacetime. A null 



hyperplane cut can also be given in terms of three de- 
grees of freedom: the null normal to the hyperplane (for 
which there are two degrees of freedom, the anti-celestial 
sphere) , and a parameter describing the translation of the 
hyperplane away from the vertex of the cone. This last 
degree of freedom also exists for spacelike hyperplanes, 
but because the intersection of the spacelike hyperplane 
with the null cone is compact (specifically a two-sphere) , 
one can fix this translation degree of freedom by fixing 
the area of the sphere. In the case of a cut by a null hy- 
perplane, the intersection is noncompact, so this degree 
of freedom cannot be fixed. 

Though the degeneracy measure 0^- may be no more 
well-defined in general than Ay- , there are still reasons to 
prefer it for purposes of defining a notion of approximate 
Petrov class. The main reason is that when a null hyper- 
plane cut is made through a null cone, one special null ray 
is singled out: the one parallel to the hyperplane. Again, 
the intersection of the future null cone with a null hy- 
perplane is itself a spacelike two-dimensional plane, and 
because the null ray parallel to the hyperplane never in- 
tersects the hyperplane, it is only represented on the in- 
tersection plane as a point at infinity. Eq. (|47|) shows that 
the null ray that gets mapped to the point at infinity is 
the one that points along the tetrad n vector. This is the 
behavior that we saw in Sec. IIII1 The quasi-Kinncrsley 
tetrad naturally adapts itself to the principal null direc- 
tions in the ringdown to Kerr geometry, such that two 
of them fall toward the origin of the C plane and two 
approach infinity. This is because the quasi-Kinnersley 
tetrad is designed to adapt itself to nearly-degenerate 
principal null directions. To the extent that the nu- 
merical tetrad approximates a quasi-Kinncrsley tetrad (a 
common implicit hope for the extraction of gravitational 
waveforms) this behavior will be seen also in numerical 
simulations. An example of this will be seen in the next 
section. 



V. NUMERICAL RESULTS 

Our numerical implementation of these mathematical 
tools begins, as in Ref. [f|, with the fourth-order polyno- 
mial in Eq. §5§ . We begin by computing the Weyl scalars 
in a reference tetrad. The timelike orthonormal tetrad leg 
e*o is taken to be the normal to the spatial slice, and the 
spacelike orthonormal tetrad legs are constructed from a 
Gram-Schmidt orthogonalization of the basis vectors of a 
spherical-like coordinate system within the spatial slice, 
essentially similar to the method in [f|. This tetrad is 
singular at the z axis, as the complex phase of the rh leg 
becomes undefined due to the coordinate singularity of 
the spherical chart, but all quantities we present will be 
independent of this complex phase, and thus will have 
well-defined values on the axis. 

Our code computes the electric and magnetic parts 
of the Weyl curvature tensor directly from data on the 
spatial slice, using Gauss-Codazzi relations and assuming 



10 



the Einstein equations are satisfied and that no matter 
fields are present: 

Eij = ( 3 Rij+KK i:i -K ik K^ STF (50) 
R — (f mn r> K -\ fRFl 

Here, Kij is the extrinsic curvature of the spatial slice, 
Di is the torsion free covariant derivative compatible with 
the spatial metric, Rij its Ricci curvature, e^-fc the spa- 
tial Levi-Civita tensor, and the subscript STF means that 
the quantity in brackets is made symmetric and tracefree 
in the indices i and j. Once these tensors are computed, 
we construct the Weyl scalars from them as in Eqs. (4) of 
Ref. [lj|, using the radial tetrad leg for u, and the polar 
leg plus i times the azimuthal leg for m. 

Once the Weyl scalars are known, one can go about 
solving for the roots Xi of Eq. ([3]). We do so point by 
point on the computational grid with simple Newton- 
Raphson iteration and polynomial deflation ^7[- In 
many cases, these methods are conventionally followed 
by root polishing — using the computed roots of the 
deflated polynomials as initial guesses in new Newton- 
Raphson iterations of the initial polynomial, with the 
hope of correcting roundoff error accumulated in the de- 
flation process — but in this case root polishing has no 
noticeable effect. This is presumably because the roots 
under consideration are very nearly degenerate, so error 
in the Newton-Raphson iterations themselves dominates 
the error accumulated in the polynomial deflation. 

As in Ref. @, we focus our attention on a simulation 
of the ringdown of a binary black hole merger to Kerr ge- 
ometry. The simplest example of such a merger is one fol- 
lowing the inspiral of two equal mass nonspinning black 
holes in a noneccentric configuration. This data set was 
presented in detail in Ref. [l8[ , and the multipolar struc- 
ture of the post- merger horizon was studied in Ref. fl9| . 
In the former paper, it was noted that two independent 
measures of black hole spin, designed to agree if the final 
black hole is Kerr, agreed to well within the estimated ac- 
curacy of the numerical truncation. In the latter paper, 
this correspondence was studied in much greater detail, 
demonstrating that all of the multipole moments on the 
apparent horizon that we were able to compute agreed 
very well with those of the Kerr horizon (see Refs. [20l.l2i| 
for a similar analyses). While these provide a very com- 
pelling case that the final black hole is Kerr, they do not 
present a major benefit of the methods described here 
and in Ref. @: being fully local. The degeneracy mea- 
sures described here can be computed independently at 
each point in spacetime, rather than simply for the ap- 
parent horizon as a whole. In this way one can imagine 
demonstrating not just the fact that a spacetime is set- 
tling down to Kerr geometry, but where it is doing so 
more quickly and more slowly, and possibly even the re- 
lationship between the approach to Kerr geometry and 
the presence of gravitational radiation. 

This locality of the approximate Petrov classification 
system, while beneficial for the reasons described above, 
unfortunately comes at the cost of another type of gauge 
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FIG. 4: The ringdown after the merger of two equal mass, 
initially nonspinning holes. The curves show the behavior of 
the two smallest values of each of the degeneracy measures 
Aij and <dij, with respect to coordinate time, evaluated at 
z = 4.5, x = y — in an asymptotically inertial coordinate 
system. The heavier curves show the two values of the <dij 
measure, the lighter curves the Aij measure. The solid curves 
show the smaller values of these measures, and the dashed 
curves the larger. The symmetries along this axis force the 
tetrad to satisfy the basic conditions of a "quasi-Kinnersley" 
tetrad, as described in section IIHI and the results of that 
section explain the initial exponential growth of the higher 
An curve. 



ambiguity. If one wishes to investigate the time depen- 
dence of the degeneracy measures, then one must choose 
a worldline in spacetime along which to compute these 
quantities. In principle one could reduce this ambiguity, 
for example by computing along timelike geodesies, or 
worldlincs preferred by some sort of symmetry, if any ex- 
ist. For example, the symmetries inherent in a merger of 
equal mass, initially nonspinning holes provide at least 
one preferred axis for consideration. The initial data 
satisfy a discrete symmetry under 180-degree rotations 
about a certain axis, taken in our simulations to be the 
coordinate z axis, along which the initial "orbital angu- 
lar momentum vector" can be intuitively said to point. 
To the extent that the numerical simulation preserves 
this discrete symmetry, the z axis sweeps out a geomet- 
rically well-defined worldsheet in spacetime. In principle 
this timelike worldsheet could be restricted to a single 
well-defined timelike worldline, on which data can be ex- 
tracted, by intersecting the worldsheet with a level sur- 
face of some curvature invariant. Here, however, we do 
not go to such lengths, electing instead to follow coordi- 
nate worldlines, as in Ref. |(|, but paying special atten- 
tion to the symmetry axis. 

In Fig. [H data are shown for the two smallest — small- 
est among the various possible root pairings — values of 
the degeneracy measures Ay and O^ evaluated at the 
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coordinate location x = y = 0, z — 4.5 7 as a function 
of coordinate time after the formation of the common 
apparent horizon in the datasct described in Ref. fl8| . 
Because the tetrad is, by construction, adapted to the 
symmetry axis at this location, 8 it is forced to be "trans- 
verse" in the sense of Sec. IIIII (a fact which we have ver- 
ified by a direct inspection of the computed values of 
\^i\ and | ^3 1) . In Fig. @J we initially see exponential 
growth in the second-smallest root difference Ay , as one 
would expect from the considerations of Sec. IIIII Even- 
tually, this exponential growth gives way to exponential 
decay, similar to the behavior seen in Fig. [3] This occurs 
because the data we compute here are actually interpo- 
lated to the polar axis from data on grid points slightly 
offset from it. On these grid points, the tetrad differs 
slightly from the quasi-Kinncrsley tetrad, as in the dis- 
cussion near the end of Sec. IIIII One might hope that 
this eventual decay would only occur on these offset grid 
points, and that the data interpolated to the axis would 
grow indefinitely as \/— 6^2/^4, but as the black hole 
settles down the growing peak in Fig. [5] shrinks in width, 
so eventually one would expect it not to be resolved by 
the spectral discretization. Incidentally, we have con- 
firmed that the rates of exponential decay in the decaying 
curves, and the rate of exponential growth in the grow- 
ing curve, each roughly equal half of the damping rate 
of the least-damped quasinormal mode of a Kerr black 
hole of the same final mass and spin as our final rem- 
nant. One would expect this from Eq. (|29p . The most 
important issue to note about Fig. [U though, is the dis- 
crepancy between the picture implied by the Ay values, 
and that implied by the Qij values. The highest and low- 
est curves are the two relevant values of Ay. At early 
times, even the qualitative behavior of these curves are 
different, one growing and one decaying. Even at late 
times, when both curves decay exponentially (and even- 
tually settle to fixed limits due to numerical truncation 
error), they still differ by four orders of magnitude. If 
Ay were naively interpreted as defining the "nearness" 
to any Petrov class, then the response to Fig. [4] would 
be that the final result of the numerical simulation is of 
Petrov Type II, not Type D, on this axis. The other 
two curves in Fig. [4] tell a very different story. The two 



For a sense of scale we note that the apparent horizon, in these 
coordinates, settles down at late times roughly to a coordinate 
sphere with radius of 2.61. 

Actually the tetrad is, strictly speaking, not well defined on the 
axis, because it is constructed from the spherical coordinate ba- 
sis, which is singular there. However, for the objects we com- 
pute, this singularity has no effect. The £ and n tetrad legs arc 
well-defined on the axis, it is just the complex phasing of the fa 
vector that becomes undefined there. The actual quantities we 
compute, however, Ajj and &ij, are invariant under spin trans- 
formations (spin-boost transformations with |c| = 1 in Eq. 1151 1 ). 
Because there are no grid points on the axis, these quantities can 
always be computed, and because they are spin invariant, they 
can be smoothly interpolated to the axis. 
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FIG. 5: The two smallest values of the degeneracy measures 
Aij and Qij evaluated at the coordinate location x = 4.5, y = 
z — in the ringdown after a merger of equal-mass, initially 
nonspinning black holes, as in Fig. [4] The lighter curves, 
which are generally the highest and lowest curves, are the 
two values of the Aij measure. The heavier curves represent 
the Qij measure. Initial growth of the larger A^ value is still 
present, but less striking than in Fig. [4] Nonetheless, the two 
A^ values again differ by multiple orders of magnitude, while 
the two Qij values generally differ by only one. 

heavier curves in Fig. |4] show the two smallest values of 
the Qij measure. Both curves fall exponentially at the 
same rate, and they lie within roughly a factor of ten of 
one another throughout the entire ringdown. According 
to Qij the spacetimc falls quite unambiguously to Petrov 
Type D on the polar axis. 

The behavior at different coordinate locations is less 
striking, but shows roughly similar features. Fig. [5] 
presents the same quantities as Fig. @] computed instead 
at x = 4.5, y = z = 0. Again, the highest and lowest 
curves are the two relevant values of Ay . The higher one 
grows slightly (on average) in the early ringdown, but 
settles into exponential decay quite a bit sooner than in 
Fig. [31 and throughout the ringdown remains separated 
from the values in the lowest curve by roughly two to 
three orders of magnitude. This still, however, provides 
a marked contrast from the two 0y curves, which again 
lie within roughly a single order of magnitude of one an- 
other throughout the ringdown. A particularly interest- 
ing feature is visible in the uppermost curve of Fig. [5] from 
the beginning of the post-merger dataset to roughly the 
coordinate time 8275. This time frame is magnified in 
Fig. [6] Because the spacetime is symmetric under reflec- 
tions across the z = plane, the spatial projections of 
the principal null directions on this plane must either be 
tangent to the plane, or otherwise reflection-symmetric 
across it. The jagged peaks in Fig. |6] imply that the for- 
mer possibility seems to apply here. When the tetrad 
ft vector happens to point along a principal null direc- 
tion, the corresponding C value of the polynomial root 
becomes infinite. If the spatial projections of two of the 
principal null directions lie in the same plane as the spa- 
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FIG. 6: Magnification of the highest curve in Fig. [4] The 
sharp peaks in the bottom curve before time 8275 imply that 
principal null directions are occasionally crossing the tetrad 
legs. The upper curve corroborates this by explicitly showing 
that the absolute value of the Weyl scalar $4 approaches zero 
at times coincident with these sharp peaks. 



10'r 




coordinate time (code units) 



FIG. 7: L2 norms of the two smallest values of Aij and 
Oij, integrated over a thick spherical shell extending from 
just within the apparent horizon, r — 2.2 in code coordinates, 
to an outer boundary at r = 5 code coordinates. Again, 
both values of <dij fall quickly to zero at the same rate, while 
the larger Ay value hangs up initially, and eventually falls 
only to a level over two orders of magnitude above its smaller 
counterpart. 



tial projection of the tetrad n vector, and they oscillate in 
direction, repeatedly crossing the spatial projection of n 
(due either to physical or gauge effects), then one would 
expect the A. y - values involving these principal null di- 
rections to show sharp, repeating peaks, such as those 
in Fig. [SI Eventually, such crossings could be expected 
to stop as the angle between the spatial projections of 
the principal null directions becomes small and as gauge 
and physical oscillations become less wild. After roughly 
t = 8275 in these code units 9 the oscillations in this curve 
become more smooth, implying that the principal null di- 
rections are no longer crossing n. 

Figure [7J presents L2 norms of the same degeneracy 
measures for the same ringdown dataset. While this 
avoids the choice of a specific coordinate location for 
observation, we should note that there is still a certain 
amount of coordinate ambiguity in this quantity. The 
L2 norm is only over a certain subset of the spatial slice. 
The inner boundary is the excision boundary of the sim- 
ulation, slightly inside the apparent horizon. The outer 
boundary is a boundary of subdomains in the simula- 
tion, with coordinate radius 5. The purpose of this outer 
boundary is to avoid numerical difficulties when the Weyl 
scalars become too small to calculate accurately the roots 
of the polynomial in Eq. ©. Again in Fig. [7J we see the 
larger value of Ay hanging up while both values of 9, 3 - 
decay exponentially. 



To aid in translating the code units to physically relevant units, 
we note that the final mass of remnant black hole, in these code 
units, is roughly 1.98. 



VI. DISCUSSION 

The primary goal of this paper has been to explain the 
peculiar behavior, noted in Ref. @, that the spacetimes 
of binary black hole mergers seem to "hang up" in Petrov 
Type II, and if they fall to Type D at all (according to 
one's choice of tolerance), they do so much later. 

Properly clarifying this issue has required us to in- 
vestigate the geometrical meaning of the polynomial in 
Eq. ^ and the measure of degeneracy principally used 
in Ref. [6[ , the absolute value of the difference of any two 
roots, := |A, — The true space of interest is the 
space of future-directed null rays at a point, the so-called 
anti-celestial sphere. As argued in Sec. IIV1 the complex 
quantity A acts as a coordinate on this two-dimensional 
space. It should therefore not be surprising that A^, a 
coordinate distance, fails to represent geometries in the 
space of null rays, since it is impossible to cover a topo- 
logical sphere with a single coordinate chart without the 
presence of coordinate singularities. 

It would therefore seem that the right thing to do 
would be to consider truly geometrical distances in the 
space of null rays as defining the nearness of principal 
null directions to one another. This approach, unfor- 
tunately, is clouded by the nonexistence of a preferred 
metric structure on the anti-celestial sphere. The anti- 
celestial sphere has a six-dimensional conformal group, 
corresponding to the proper Lorentz group of Minkowski 
spacetime. While this group carries a three-dimensional 
subgroup of isometries — corresponding to rotations — 
which have no effect on the "distances" between any two 
null rays, the three remaining dimensions — correspond- 
ing to boosts — conformally rescale the metric on the 
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space of null rays. For this reason, fixing a unique ge- 
ometry on the space of null rays requires fixing a unique 
"observer" with respect to which this boost freedom is 
fixed. In the introduction, we described similar difficul- 
ties in attempting to define a concept of approximate 
Petrov class by algebraic means. 

In Sec. II VI we also presented a geometrical interpreta- 
tion of the quantity Ay. Rather than simply as a co- 
ordinate distance on the space of null rays, Ay can be 
interpreted as a metric distance along the- cut of a null 
cone made by a null hyperplane. In a sense, one is here 
adapting the metric on the space of null rays to a null 
observer. Similarly, is a distance function on the in- 
tersection of the future null cone with a spacelike plane 
orthogonal to our timelike observer. 

There are equally many degrees of freedom in cutting 
the null cone with a spacclikc hyperplane as there are in 
cutting it with a null hyperplane (once one restricts the 
spacelike cuts to normalize the area of the anti-celestial 
sphere, a restriction that cannot be made on null hyper- 
planes because the intersection is noncompact). For this 
reason the degeneracy measure that we have introduced, 
&ij in Eq. (|45|) , cannot be considered fundamentally any 
more invariant than Ay, though in practice it is easier 
to imagine a fleet of preferred timelike observers than of 
null observers, such as stationary observers in stationary 
spacetimes, observers adapted to timelike approximate 
Killing vectors in approximately stationary spacetimes (if 
such vectors can be reasonably defined), or freely falling 
observers following timelike geodesies from infinity. We 
have not, however, attempted to find any such preferred 
classes of observers in the numerical results presented 
here, either for fixing the geometry on the space of null 
rays or for fixing the worldlines along which data are 
extracted. The main value of the new degeneracy mea- 
sure 9y is not that it is more gauge-invariant, but rather 
that it naturally captures the compactness of the space 
of null rays, and thereby avoids relegating any particular 
null ray to a point at infinity. 

The need to avoid relegating any null ray to a point 
at infinity is particularly acute in practice, as the rays 
at infinity in the physically preferable quasi-Kinnersley 
tetrads become the principal null directions themselves 
as a spacetime settles down to Kerr geometry. This be- 
havior was investigated in Sec. IIIII As principal null di- 
rections relax to degeneracy at the point at infinity in A 
space, the degeneracy measure Ay grows exponentially 
rather than decaying exponentially as one would naively 
expect. While the tetrads used in numerical simulations 
are not commonly quasi-Kinnersley tetrads, there is gen- 
erally an implicit hope, for purposes of wave extraction, 
that they are close to it in some rough sense. Indeed, 
as is visible in Fig. SI this nonintuitive behavior in the 
quasi-Kinnersley tetrad can corrupt measurements of Ay 
in even a simple coordinate-adapted tetrad (c/. Fig- [3j) . 



The other figures in Sec. |V] tell much the same story, 
though in somewhat less dramatic terms. Fig. [5] shows 
indications of principal null directions directly crossing 
the tetrad null vectors, repeatedly causing the null direc- 
tions to be represented by the point at infinity in A space, 
due purely to the choice of spatial tetrad. Fig. [7] shows 
that the hangups in the degeneracy measure Ay are not 
limited to the partially symmetry-preferred x and z axes. 
In fact, Ay, computed everywhere in a tetrad adapted 
to spherical coordinates, clearly shows this hangup in 
Petrov Type II even in an integral L2 norm, while 0y 
does not. 

Another motivation of this paper has been to provide 
further evidence that the final remnants of our black hole 
merger simulations are Kerr black holes. This was indeed 
the central focus of Ref. 0, and they even went so far 
as to demonstrate that their final black hole has no NUT 
charge or acceleration (as in the C-metrics; see Ref. [H|), 
an issue that we have not explored here. 

The question of whether a black hole is "settling down 
to Kerr" can be attacked at a variety of levels. In a re- 
cent paper [l{|, we studied the question at a quasilocal 
level, verifying that the quasilocal source multipole mo- 
ments of the black hole settle to the values required by 
the Kerr geometry (see also Refs. [10, [H|). More recently, 
Ref. [23| presented theoretical tools for approaching the 
question at a global level (global on any given spacelike 
slice). The approach taken in Ref. @ was in part local 
(measurement of Ay), and in part global. The method 
by which Campanclli et al. verified the vanishing of the 
NUT charge and acceleration involved limits of curva- 
ture invariants to large radii. If one wishes to rule out 
NUT charge or acceleration at a local level, to provide 
a more cohesive picture when combined with local alge- 
braic degeneracy measures, this can be done with quan- 
tities presented in Ref. [ll|, though as described in the 
introduction of this paper, collapsing tensorial quanti- 
ties to scalar quantities would require a positive-definite 
background metric, which could require fixing a slicing 
or a threading of spacetime. 



Acknowledgments 

I thank Mark Schccl for providing the numerical sim- 
ulation data studied in Sec. [V] I also thank the Cal- 
tech and Cornell numerical relativity groups, particu- 
larly Jeandrew Brink, Larry Kidder, Geoffrey Lovelace, 
and Saul Teukolsky for frequent discussions. This 
work is supported in part by grants from the Sherman 
Fairchild Foundation to Cornell, and by NSF grants 
PHY-0652952, DMS-0553677, PHY-0652929, and NASA 
grant NNX09AF96G. 



14 



[1] A. Z. Petrov, Gen. Relativ. Gravit. 32, 1665 (2000). 

[2] F. A. E. Pirani, Phys. Rev. 105, 1089 (1957). 

[3] H. Stephani, D. Kramer, M. Maccallum, C. Hoense- 
laers, and E. Herlt, Exact Solutions of Einstein's Field 
Equations (Cambridge University Press, Cambridge, UK, 
2003), ISBN 0-521-46136-7. 

[4] L. Bel, Cah. de Phys. 16, 59 (1962). 

[5] R. Penrose, Ann. Phys. 10, 171 (1960). 

[6] M. Campanelli, C. O. Lousto, and Y. Zlochower, Phys. 
Rev. D 79, 08401 2 (2009), arXiv:0811.3006. 

[7] SpEC, URL |http : //www . black- holesTo rg/SpEC . html | 

[8] R. Penrose and W. Rindler, Spinors and Space-Time: 
Vol. 1, Two-Spinor Calculus and Relativistic Fields 
(Cambridge University Press, Cambridge, UK, 1987), 
ISBN 0521337070. 

[9] R. Penrose and W. Rindler, Spinors and Space-Time: 
Vol. 2, Spinor and Twistor Methods in Space-Time Ge- 
ometry (Cambridge University Press, Cambridge, UK, 
1988), ISBN 0521347866. 
[10] J. Baker and M. Campanelli, Phys. Rev. D 62, 127501 
(2000). 

[11] J. J. Ferrando and J. A. Saez, Class. Quantum Grav. 26, 
075013 (2009), arxiv:0812.3310. 

[12] A. Nerozzi, C. Beetle, M. Bruni, L. M. 
Burko, and D. Pollney, Phys. Rev. D 
72, 024014 (pages 13) (2005), URL 



http : //link . aps . org/abstract/PRD/v72/e024014 

[13] E. Newman and R. Penrose, Journal of 
Mathematical Physics 3, 566 (1 962), URL 
|http : //link . aip . org/link/? JMP/3/566/1 | 

[14] C. Beetle and L. M. Burko, Phys. Rev. Lett. 89, 271101 
(2002), gr-qc/0210019. 

[15] J. Stewart, Advanced General Relativity (Cambridge Uni- 
versity Press, Cambridge, UK, 1993), ISBN 0521449464. 

[16] A. Nerozzi, Phys. Rev. D 75, 104002 (2007). 

[17] W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. 
Flannery, Numerical Recipes: The Art of Scientific Com- 
puting (Cambridge University Press, 2007), 3rd ed. 

[18] M. A. Scheel, M. Boyle, T. Chu, L. E. Kidder, 
K. Matthews, and H. P. Pfeiffer, Phys. Rev. D 79, 024003 
(2009), arxiv:0810.1767. 

[19] R. Owen, Phys. Rev. D 80, 084012 (2009), 
arxiv:0907.0280. 

[20] E. Schnetter, B. Krishnan, and F. Beyer, Phys. Rev. D 
74, 024028 (2006), gr-qc/0604015. 

[21] M. Jasiulek, Class. Quantum Grav. 26, 245008 (2009), 
arxiv:0906.1228. 

[22] J. B. Griffiths and J. Podolsky, Exact Space-Times 
in Einstein's General Relativity (Cambridge University 
Press, Cambridge, UK, 2009), ISBN 978-0-521-88927-8. 

[23] T. Backdahl and J. V. Kroon (2010), arxiv:1001.4492. 



